Energy distribution of maxima and minima in a one-dimensional random system 
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We study the energy distribution of maxima and minima of a simple one-dimensional disordered Hamiltonian. We find that 
in systems with short range correlated disorder there is energy separation between maxima and minima, such that at fixed 
energy only one kind of stationary points is dominant in number over the other. On the other hand, in the case of systems 
with long range correlated disorder maxima and minima are completely mixed. 
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When the statistical properties of a system are studied, 
great attention is usually devoted to its ground state and 
to the first excited states. Moreover, in the case of ran- 
dom systems it is known that, in addition to the lowest 
energy states, also metastable states are important, es- 
pecially for the dynamical evolution of the system. As 
a consequence, the physical understanding of disordered 
models is most of the time founded on the properties 
of absolute as well as local minima of the Hamiltonian, 
while the role of stationary points of different nature is 
in general disregarded. In spite of this, the conviction is 
growing up that stationary points different from minima 
do have an importance, both from a dynamical and a 
static point of view. 

Many different disordered systems display an off- 
equilibrium dynamical behavior which is suitable to be 
interpreted in terms of non-trivial structure of their phase 
space. Among these we find structural glasses [jlj, spin 
glasses [^|^], random manifolds Q and neural networks 
In all these cases the geometric structure of the en- 
ergy landscape is often invoked in order to give at least 
an intuitive picture of the relaxational dynamics. In this 
context it is clear that also unstable stationary points 
must be taken into consideration. For instance, the pre- 
sence of fiat directions in the phase space, marking a 
borderline between the last stable minima and the first 
unstable saddles, has been proposed as a possible expla- 
nation of slowness in glassy systems [|]-|| . 

In the light of these considerations, we believe it is im- 
portant to understand the connections between the phys- 
ical properties of a disordered system and the geometri- 
cal structure of all the stationary points of its Hamil- 
tonian. Unfortunately, it is in general very difficult in 
an TV-dimensional model to discriminate the stationary 
points according to their degree of instability. In this Let- 
ter we will thus focus on a simple one-dimensional case 
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and exactly compute the average energy distributions of 
maxima and minima. We will find a simple connection 
between the nature of the disorder which rules the physics 
of the system and the mutual distributions of the station- 
ary points. Despite its simplicity, we expect the model 
studied here to capture at least some of the main features 
of more general problems. 

Let us consider the one-dimensional random Hamilto- 
man §0, 

H(x) = ^mx 2 + V(x) , 

where the position x is a real variable and the mass to is 
a parameter. V(x) is a Gaussian random potential, with 
zero average and variance V{x\)V{x2) = G{x\ — x 2 ), with 
G(x) = G(—x). The statics and the dynamics of this 
model have been studied both for the one-dimensional 
case jjTcj and for the more general iV-dimensional case 

IMF 

The number of stationary points of H is determined 
by the competition between the random potential and 
the harmonic mass term. This number is large for 
small to, whereas only one single minimum is present 
at large m. The physical properties of this model are 
encoded in the function G. In order to understand its 
meaning we consider the average displacement [A(d)] 2 = 

[V{xi) - V{x 2 )f = 2G(0) - 2G{d), where d = (x± - x 2 ) 
is the distance. Once introduced A it is natural to define 
two different classes of random potentials. If A(d) goes 
to a finite value A(oo) for d — > 00, then the memory is 
lost after a finite distance and V is called short range 
(SR). On the other hand, if A(d) ~ <P (7 > 0), then the 
displacement grows indefinitely with d and the potential 
is long range (LR). In the SR case we can assume with- 
out loss of generality that G(x) is a positive even function 
which is zero at infinity, so that A (00) = y/2G(0). In the 
LR case we have to be more careful, since a diverging dis- 
placement would require G(d) — > —00 for d — > 00, which 
is incompatible with the condition of having a positive 
kernel in the functional distribution of V. In order to 
correctly define the LR model we must put the system 
in a box of size L and define Gr,(x) through its Fourier 
transform, 
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G L {x) = - [ dqG(q)e^ x 
T Jl/L 



(1) 



The function G{q) must be positive and for the LR case 
must be not integrable in zero. In order to avoid any ul- 
traviolet divergence we can assume both for the SR and 
LR cases G(q) to decay at infinity faster than any power. 
We can thus define these two classes of models simply in 
terms of the behavior of G at zero momentum. In this 
way the LR case is well defined: Al(g?) increases indef- 
initely with d and Vl{x) 2 — Gl(0) diverges with L, as 
expected, since in a LR random potential the uncertainty 
on the height V of one single point x increases with the 
size L of the system, while it remains finite in a SR po- 
tential. Our analysis will not depend on the explicit form 
of G. 

Let us denote by Af/f (E, to) dE, the number of statio- 
nary points of H(x) with degree of instability equal to 
k (k — for minima, k — 1 for maxima), which have 
energy between E and E + dE, for a given mass m. The 
superscript V indicates that this distribution corresponds 
to the sample V. Eventually we shall average over V. 
The distribution A/% (E,m) is given by, 

N%(E,m)= J dxS(H')S(H - E)\H"\5[9{-H") - k] . 

In order to handle the modulus and the ^-function we 
use the following relations, 

0(~H") = -L km [\og(H" - ie) - \og{H" + ie)] , 
\H"\ = lim(H" + ie)?(H" -ie)i . (2) 



Using an integral representation for the (5-function, we 
can write, 

1= \H"\5[6{-H") -k] = 

J dfie 1 ^ 1 (H" + ie)?+& (H"-ie)*-£. 

The last two factors can be rewritten using the identity, 
(H"±ie) n ± = J dx b ± d X b ± exp ^-f^ ± {H ff ±ie)xi\ , 

where x± an< ^ X± are Grassmann variables and the an- 
alytic continuation n± — > (1/2 ± /J,/2ir) must be done. 
As a next step we define the Grassmann vector [ p"5[ : 
if> a = (x\ ■ ■ ■ X+ + , X — * * ■ X™~ ), which allows us to write, 

1 = J dn e lk ^ J #„ # a exp ^- J2 M H " + , 

where the vector e a is split into two parts: e a = e for 
a < n+, e a = —£ for a > n+, and n = (n+ + n~) — > 1. 



Note that this replica approach can be easily generalized 
to N dimensions. 

Let us introduce in the expression for A/% the Lagrange 
multipliers A and w, to represent respectively d(H') 
and 5(H — E) . The F-dependent part then becomes: 
exp[(iu> + iXd x — ip a iJ} a d x d x )V(x)], which can be averaged 
over the Gaussian distribution of V. This produces a 
quartic term (X) a ^a^a) > that can ^ e ma de quadratic 
by means of a Hubbard-Stratonovich transformation, in- 
troducing an auxiliary variable y. It is now possible to 
perform all the Gaussian integrals over (X,x,ip). This 
gives a term (m + y + ie) 1 / 2+/i / 27r (m + y — ie) 1 / 2 ^^/ 2 ^ , 
which, using back relations can be written as \m + 
y\ exp[—ifj,9(—m — y)]. Integrating over /i we finally ob- 
tain the average distributions Afk(E, m) = A/% {E, to), 

/+oo 
dyF(y,E,m), (3) 
-m 

/ — m 
dyF(y,E,m), (4) 
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where a = (7(0), ai = -G"(0) and a 2 = G""(0). The 
difference T>(E,m) — Afo (E,m) — Af\ (E, to) between mi- 
nima and maxima has a much simpler expression, 



„- , t-, s 1 f duj , 

V(E, m) = — / — VroT 



By integrating equations (y) and (|j) over the energy we 
get the total number of minima and maxima, A/o(m) 
and A/i(m), at a given value of the mass. Note that, 
as required by the Morse theorem JT^j] , the total num- 
ber of minima minus maxima is equal to one, that is 
Af (m) -A/i(m) = J dE V{E,m) = 1. The explicit ex- 
pression for the total number of minima is, 
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Afo(m) is a smooth function of to which goes to one when 
to goes to infinity and starts increasing very steeply at 
masses smaller than m ~ y/av- This value of the mass 
marks a crossover from the region where only one min- 
imum exists to the region where many different minima 
(and maxima) appear. As expected this mass is the 
same critical mass as the iV-dimensional mean-field case, 
iTLc = \J 0,2/3, where a glassy transition occurs [ ^0|Jl~2]| . In 
the following we will always consider to < m c . 

We analyze now our results, starting with the SR 
potential. In Fig.|l| we plot Afo(E,m), Afi{E,m) and 
T>{E, to) functions of the energy E, for m < m c . 
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FIG. 1. SR potential: M Q {E,m), Afi(E,m) and T>(E,m), 
functions of the energy E, at m < m c . Here G(q) = e~ q . 

The first thing we notice is that the two curves of mi- 
nima and maxima are quite separated one with respect 
to the other, so that their peaks do not overlap. As a 
consequence T>(E, m) gives by itself a rather clear pic- 
ture of the distribution of the different stationary points 
and at low energies it approximates well Mq(E, m). This 
is important because T> is always very simple to com- 
pute, being V(E,m) = J dxS(H') H" S(H - E). Thus, 
the computation of this quantity does not require the 
modulus, nor the ^-function, which are in general very 
difficult to treat. In other words, in the SR case there is 
a partial decoupling between maxima and minima, which 
is sharper the lower is the energy. As a consequence, at 
fixed energy only one kind of stationary points is dom- 
inant over the other and No(E,m) ~ T>(E,m), for low 
enough E. It is remarkable that this holds for the SR 
potential. Indeed, it has been proved in Q that in the 
Af-dimensional p-spin spherical spin glass, which belongs 
to the SR class p^] , an identical phenomenon occurs: in 
that mean-field model, besides minima and maxima there 
are saddles of any order, but at fixed energy only one 
kind of stationary points is dominant over the other, so 
that the number of minima of the mean-field free energy 
(and therefore of states) can be safely calculated via the 
approximation Afo(E) ~ T>{E). We note that this same 
approximation has been used many times in the context 
of spin glasses, regardless of its grounding p6| . 

Let us now turn to the LR potential. We stress that 
the model is defined as long range (and thus it is different 
from the SR case) only in the limit L — > oo. From Eq.(Q) 
we know that ao = Gl(0) = V(x) 2 diverges with L. As 
mentioned above, the physical meaning of this is that 
the uncertainty on the value of H in any point x diverges 
when L goes to infinity. As a consequence, the energy E 
is no longer a good variable to label the height of the sta- 
tionary points. In order to keep everything well defined 
in the limit L — > oo it is thus necessary to measure the 



energy in units of the natural diverging scale, that is i/ag. 
Therefore, we must define a rescaled energy £ = E/^/ao 
and study the distributions of maxima and minima as 
functions of £ . Denoting these new distributions by Po 
and Pi, we have, Pk(£, m) d£ = Nk{E, m) dE. We stress 
that £ is the only variable we can sensibly regard as the 
energy for the LR potential. Note, on the other hand, 
that this rescaling is irrelevant for the SR case, where ao 
is finite. Taking the limit L — > oo in Eq.(|^) we find, 

P .i(£,m) -> A/o.i(to) -4= e"^ 2 . (6) 

V 27T 

This equation shows that in the LR case the two distri- 
butions are just the same function, scaled by the total 
number Ao(m) or Miirn). Maxima and minima are no 
longer separated in energy. Indeed, for m <C m c , we have 
A/o (to)/A/i (m) ~ 1 and the two curves collapse one on to 
the other. The conclusion is that when the total number 
of stationary points is large in a LR system, maxima and 
minima are completely mixed together, so that at each 
given energy they are equally numerous. Thus, in stark 
contrast with the SR case, no decoupling of the station- 
ary points occurs, no matter how low is the energy. 

A further step is necessary to prove that this mixing 
in the LR case is a typical behavior and not simply an 
artifact coming from the average. Indeed, it is possible 
to think of a system where sample by sample maxima 
and minima are well separated, but where the mixing 
described above appears only after averaging over differ- 
ent samples. As an example, we consider the family of 
Hamiltonians H w (x) = sin(x) + w, where w is a random 
variable with zero average and variance a. It is clear that 
for each sample maxima and minima are perfectly sepa- 
rated, since Mq 1 (E) = S(E±1 — w). However, averaging- 
over w we obtain two distributions with separation be- 
tween their averages equal to 2 and variance a. Thus, if 
we rescale the energy by a factor y/a and take the limit 
a — > oo, we would conclude that there is mixing between 
maxima and minima, which is sample by sample false. 

In order to prove that the LR potential does not cor- 
respond to such an artifact, we consider the statistics 
of the extreme values of the Hamiltonian. Let us de- 
fine the two distributions Aq(E) = S(E — Emin) and 
Ai(E) = 5(E - E M ax), where Emin and Em ax are 
the energies of the absolute minimum and maximum of 
H(x) (we consider as absolute maximum the highest local 
maximum). The separation between these two distribu- 
tions is t\A = E M ax - Emin = {E)ai - (E)a , and let 
S be their variance. Consider now the ratio AA/y/S. It 
is easy to see that in the artificial case described above 
this ratio goes to zero when the variance of the disorder 
a goes to infinity, since AA is finite, whereas S diverges 
as cr. On the other hand, for the LR potential the ratio 
AA/yS remains finite in the limit L — > oo. This is sim- 
ple to prove exactly for 7 = 1/2, which corresponds to 



3 



a Brownian random potential of size L || . For the gen- 
eral case the idea is that in a LR potential both AA and 
VS diverge as U 1 @, 

as can also be checked by means 
of numerical simulations. The divergence of A A implies 
that the variance of the energy distribution of maxima 
and minima is sample by sample diverging as well with 
L. Therefore, unlike what happens in the case of the 
artifact, the average scenario of the long range poten- 
tial, where maxima and minima are completely mixed in 
energy, is the typical one. 

In this Letter we have shown that the energy distribu- 
tion of maxima and minima in a one-dimensional random 
system is radically different whether the disorder is long 
range or short range correlated. An important issue is the 
extension of this investigation to N dimensions, where we 
expect to find a qualitative similar behavior, at least at 
the mean-field level. 

Indeed, as mentioned above, the p-spin spherical spin 
glass is a clear example of an iV-dimensional mean-field 
SR system where at sufficiently low energies a decoupling 
between stationary points of different nature does occur. 
Moreover, a crucial feature of this SR model is that the 
asymptotic dynamical energy reached by the system is 
exactly the same energy where the stationary points de- 
couple 0, and is larger than the equilibrium one, so that 
the dynamics never gets to the equilibrium landscape. 

On the other hand, let us assume that for a LR TV- 
dimensional model a generalization of equation (j^) is 
valid, so that all the stationary points collapse on one 
single distribution, since saddles of any degree of insta- 
bility will be bounded in energy by maxima and minima. 
This allows us to put forward the following hypothesis: 
if, as indicated by our results, in a LR system there is no 
separation at all between different stationary points, then 
no decoupling energy can exist, which means no energy 
level capable of trapping the system. In such a situation 
we would expect the dynamics to reach the minimum 
available energy, that is the equilibrium energy. 

An evidence of this conjecture can be found in the con- 
text of mean-field models for spin glasses. Here two very 
different classes of systems exist, a first class where the 
dynamical energy is larger than the equilibrium one and 
a second class where these two energies are the same. 
What has been noted in , is that the first class corre- 
sponds to systems with SR disorder (as the p-spin model), 
whereas to the second class belong LR disordered sys- 
tems. This correspondence finds its natural explanation 
in the framework we have depicted above: the inability of 
a SR system to dynamically reach its equilibrium energy 
is due to the existence of an energy level below which 
stationary points of different nature are decoupled, while 
this cannot happen in the LR case. As we have tried to 
show in this Letter, whether this decoupling occurs or 
not, is an information encoded in the energy distribution 
of all the stationary points. 
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